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I. INTRODUCTION 



Several possible approaches exist to embed u new phenomena!'' within 11 old paradigms" . 
For instance, theoretical principles that apply to the old paradigm, or that are of great 
foundational significance, can be extended and used as guidelines to include new observed 
facts. A prime example of such principles are symmetries belonging to the old paradigm, 
or natural extensions of them. In some cases, this approach leads to including several 
unobserved components to the original theory, and hence redundancy in what is actually 
needed to interpret new observations. Another perfectly reasonable pathway is, instead, to 
pursue the idea of " 'minimality'' : what is the most economical extension to the pre-existing 
framework that allows to explain the new phenomena? While quantitatively defining the idea 
of a "minimal" extension is non-trivial, the particle content, or the number of needed new 
additional parameters are natural choices to quantify whether an extension to an elementary 
particle theory is or not economical. 

One of the most compelling reasons to explore what might lie beyond the Standard 
Model (SM) of particle physics is the mysterious nature of the dark matter that dominates 
the matter content of the universe. In this context, theories that provide a dark matter 
candidate are widely considered more interesting extensions to the SM than those models 
that fall short of providing one. Supersymmetry l| and theories with universal extra- 
dimensions Q are examples of extensions of an old paradigm (here, the SM of particle 
physics) guided by symmetry principles that permit to explain new phenomena (the particle 
nature of the dark matter). The other pathway mentioned above, instead, has also been 
pursued successfully, postulating ad hoc, minimal extensions to the SM that encompass a 
dark matter particle candidate (for a systematic approach see e.g. Ref. 3)]). 

In many respects, what is widely considered to be the simplest, if not the most economical 
choice to embed a particle dark matter candidate into the framework of the SM, is to add 
a gauge-singlet real scalar field S with renormalizable interactions only, and enforcing the 
Z2 symmetry S — > —S. As we shall detail below, this theory, assuming S is the only new 
degree of freedom at the electro-weak scale, only adds three new parameters to the SM: a 
mass term, a self-interaction term and a parameter that controls the coupling of the singlet 
to the other SM fields. 

The theory we consider here has a quite long history, having been first envisioned by 
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Veltman and Yndurain (4], who introduced a scalar U U particle" to the SM and studied the 
impact of such particle on one-loop SM radiative corrections, in particular to WW scattering. 
The theory was considered in a cosmological setup, and the scalar particle - there dubbed 
" scalar phantom" - as a dark matter candidate by Silveira and Zee in Ref. [sj. Most of the 
associated relevant phenomenology was worked out in Ref. j^J, including the calculation of 
the relic particle abundance from thermal freeze-out in the early universe, the scattering 
rate of the scalar particle off of baryons (direct detection), the effect on the SM Higgs 
decay and even the impact on the galactic cosmic-ray flux. Following that seminal work, a 
number of refined studies have considered the same, simple extension to the SM. In Ref. js] 
the general case of an arbitrary number of complex singlet scalars was considered, with an 
emphasis again on cosmology in the early universe and direct detection. The specific case 
of one real singlet scalar was examined in great detail in Ref. 7|, including collider searches 
via anomalous Higgs decay patterns, dark matter self-interactions and constraints from the 
singlet potential. 

The real scalar singlet extension to the SM was promoted in Ref. js] to the status of "New 
Minimal" SM. With the advent of the Large Hadron Collider, several studies addressed the 
phenomenology of this paradigm with colliders, including e.g. Ref. js-ll]. A real scalar 
singlet also provides the possibility that the electro-weak phase transition be strongly first 
order, as needed to produce the observed baryon asymmetry in the context of electro-weak 
baryogenesis Q, Q. In Ref. Q and [15] TeV-scale scalar singlet extensions to the SM 



were shown topotentially have important implications for the recently observed cosmic-ray 
anomalies 16l-ll8|. 

Recent exciting results from direct dark matter experiments have triggered a renewed 



interest Jn real singlet scalar dark matter, that was invoked to interpret the DAM A 
CDMS 



20 



19|, 



24j and other anomalous signals like those detected with CoGeNT[25[. The use 



of gamma-ray data, especially from the recently and successfully deployed Fermi Large Area 
Telescope (LAT), to detect a signature from singlet scalar dark matter has also been recently 
addressed in Ref. H|- in the present study, we consider a channel that has not, to our 
knowledge, been addressed yet in this context: the pair annihilation into two, monochromatic 



high-energy gamma rays (see, however, Ref. 



28j). This channel is particularly relevant, given 



its unique spectral structure. This is unlike the previously considered continuum gamma-ray 
emission, a signal that could be confused with astrophysical backgrounds from e.g. emission 
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from galactic cosmic rays or from milli-second pulsars. Also, the Fermi-LAT Collaboration 
recently assessed the observational constraints on searches for this peculiar spectral feature 
in Ref. [29]. 



We find here that a particularly large and interesting region of the real singlet scalar 
dark matter parameter space for the two-photons annihilation channel is constrained by 
consideration of vacuum stability of the one-loop scalar potential, as recently studied in 



Ref. 



301 ] . We reconsider here those constraints, in light of the possibility that even though the 



electro-weak minimum might be meta-stable, its lifetime for tunneling to the true minimum 
of the one-loop potential might be much longer than the age of the universe. If this is the case, 
the stability of the electro-weak vacuum is still valid on cosmological scales. We show in this 
study that constraints form the two-photon annihilation channel are particularly relevant 
in these regions of meta-stable electro-weak vacuum, that also overlap in some cases with 
parameter space portions compatible with the tentative positive direct dark matter detection 
signals reported by CDMS. 

The outline of the present study is as follows: We first calculate, in section [Til the relevant 
SS — > 77 cross section. We then analyze the impact of the recent Fermi-LAT observations on 
the relevant parameter space in section IIIH and assess the impact for recent direct detection 
results in section IIV1 We explore the parameter space connected to meta-stable vacua in 
section [V] Finally, we outline our conclusions and summarize our results in section IVI1 



II. THE SS -> 77 CROSS SECTION 

We consider here a theory where a real scalar singlet S is added to the particle content 
of the SM. Imposing the Z 2 symmetry S — > — S to the theory, so that the singlet is stable 



and becomes a candidate for dark matter, the following 
renormalizable terms (we follow here the notation of Ref. 



^agrangian exhausts all possible 
13): 



£ = Csu + \d»Sd»S - b jS 2 - - a 2 S 2 H^H (1) 

where £sm is the SM Lagrangian and H is the SM Higgs doublet. This model adds the 
following three parameters to the SM: b 2 ^b^a 2 - We require that, at the minimum of the 
potential, the Higgs get a non-vanishing vacuum expectation value (VEV) v = 246 GeV, 
but that the singlet do not, (S) = 0, in order to avoid making the singlet-like particle 
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FIG. 1: Schematic Feynman diagram for the pair-annihilation of two singlets into two photons. 

unstable, and thus not a viable dark matter candidate, as well as to prevent the appearance 
of problematic domain walls. In the unitary gauge, = l/y/2(h + v,0) with h real, after 
symmetry breaking, the scalar potential becomes 

V(h, S) = -£- ^ + Xvh? + h 4 + hb 2 + a 2 v 2 )S 2 + ^ 4 + a 2 vS 2 h + ^S 2 h 2 , (2) 
4A 4 2 4 2 

where /i 2 < 0, A is the quartic coupling for the Higgs, and (— n 2 /X) 1 ! 2 = v. This potential 
is bounded from below, at tree level, provided that A, 64 > 0, and A64 > a 2 for negative a 2 . 
We see that the S mass is 

m s = b 2 + a 2 v 2 . (3) 

As explained in the introduction, models of this kind have been studied in the literature, 
and constraints on the parameters have been derived, mostly from dark matter direct de- 
tection experiments. The aim of the present study is to put further constraints by studying 
the pair- annihilation channel into two photons, with E 1 = ms, and by comparing with the 



291 ] . In order to do so, we study the cross 



photon lines limits obtained with Fermi-LAT 
section for the process shown in Fig. [H 
The amplitude for the process reads 

where s is the center of mass energy squared, the total decay width T = r(m/ l ) + Tg, 
with r(m/i) the decay width of the Higgs to SM particles and Ts = 8 ^ Re a/1 — 4m| / m\ 
the decay width of the electro-weak Higgs to SS. The latter vanishes if the channel is 
kinematically forbidden, i.e. if < 2ms. The annihilation cross section is given by 

M77 = ^ \M S s^ \ 2 - (5) 
In the above Equations, \Aih~^-yy\ 2 can be obtained from the result of the one-loop calculation 
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of the width T^^ 77 of the Higgs to two photons 3l|, |32] 

1 



h— >77 



167rmfc 



M 



h— >77 1 



1024vr 3 m h M, 2 



-^-1 ViV e 2 F-\ 2 



(6) 



In the annihilation process we study (see Fig. [TJ, the Higgs is, however, the exchanged 
particle and it can be off-shell. Therefore, in order to get the correct expression for | A^h— s--yy | 2 
we need to substitute in Eq. (JfJJ) with s. Thus we have 



(7) 



where % = spin- 1/2 and spin-1 identifies the particle running in the loop, N ci is its color 
multiplicity, is the electric charge in units of e, and 



Fx/2 = _ 2t [1 + (1- t )/( t )], 
F x = 2 + 3r + 3r(2-r)/(r), 



(8) 



with r = 4m 2 /s and 



f(r) 



sin 



n 2 



hi 



Z7T 



if r > 1, 
if r < 1. 



Plugging back into Eq. (j5J) we have 



77 



'32tt 3 (s-m 2 ) 2 + r 2 m 2 



2 

1 1 i 



(9) 



(10) 



where we used Mw = l/2gv. In the remainder of the paper we will consider the singlets to 
annihilate when they are non-relativistic, so that s ~ 4m|. 



III. CONSTRAINTS FROM SS -> 77 AND FERMI-LAT OBSERVATIONS 



The Fermi-LAT Collaboration has recently searched for monochromatic 7 rays in the 
range 20-300 Gev[l3], that would be produced by dark matter particle annihilation (for a 



study of the constraints from EGRET data see [33j). We indicate the resulting limits on 



(av ) 77 in Fig. 2(a) for three different representative dark matter density profiles: Einasto [34] 
(red dots), Isothermal [35| (green dots) and Navarro- Frenk- White (NFW) [36| (blue dots) - 



we refer the Reader to Ref. [17J for details on the analysis and on the specific assumptions 
for the dark matter density profiles. Given the recent interest in adopting the singlet scalar 
model to interpret signals that might be due to low mass dark matter (ms < 30 GeV), we 
thought it worth to extrapolate the Fermi limits to lower energies as well. In order to do 
so, and to be sufficiently conservative, we fixed the photon flux $ to 5 x 10 -9 cm -2 s -1 , a 
value which is in line with those given for the lowest energies in Ref. and we used the 
fact that (crv) 71 oc m|$. The resulting limits on the cross section are shown in Fig. 2(a) for 
rris < 30 GeV as crosses, with the same color scheme for the three profiles. The different 
lines adopt different values for the mass of the electro- weak Higgs and for specific values 
of the parameters in the potential that remain unspecified after fixing mj. In the solid lines 
we set 62 = 0, thus corresponding to a Lagrangian with S 2 \H\ 2 interactions plus a mass term 
for S only, while for the dashed line we fix a 2 = 0.05. Finally, the dot-dashed line shows 
the cross section corresponding to m^ = 150 GeV and the third parameter in the potential 
set to fulfill the requirement of an S thermal relic abundance equal to the universal dark 
matter density. 



It is clear from Fig. 2(a) that the cross section of Eq. (II 01) exceeds the Fermi limits only 



in a small region around the resonance, which happens for mg ~ 1/2%. This region is 



outlined with greater accuracy in Fig. 2(b) , where we plot A = (ms — m/ l /2)/m/ l versus m/j, 
after setting 6 2 = 0. The shape of the Figure can be understood as a combination of two 
factors: 

1. the non-trivial structure of the Fermi limits as a function of energy (i.e., here, as a 
function of the singlet mass) which depends on the astrophysical background and on 
the instrumental performance (e.g. point spread function and energy resolution); 

2. the fact that, with b 2 = 0, we have a 2 = m s /v 2 and the annihilation cross section 
is proportional to ,. 2 _ 2 2 . This second factor, in particular, explains the 
asymmetry of the plot. 

We point out that as we require S to have a relic density consistent with WMAP, then 
a 2 is determined as a function of ms and rrih, and it becomes very small when ms = l/2m/, 
(see e.g. [?], H])- This has the effect of canceling the resonance in our cross section, and, as 



a result we never exceed the Fermi constraints, as the dashed-dotted line in Fig. 2(a) shows. 
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We note that important constraints on the model under consideration here also stem 
from the continuum gamma-ray emission from SS annihilation. In particular, one of the 
most stringent constraints comes from Fermi-LAT observations of local dwarf spheroidal 



galaxies 



371 ] . Although a systematic comparison with the present constraints from the 



monochromatic gamma-ray emission is beyond the scopes of the present analysis, it is useful 
to compare the two constraints in a few indicative cases. Let us consider, for instance, 
the green solid line (rrih = 180 GeV, 62 = 0), and mg = 80 GeV, a model which is right 
at the level of the 77 constraint. The corresponding total pair-annihilation cross section 
can be simply read off dividing by the Higgs decay branching fraction into two photons, 
which is around 10~ 3 for = 2 x mj = 160 GeV. We thus find that (crv) tot ~ few x 
10 -25 cm 3 s _1 . This is right around what found in Ref. [stJ for a dark matter mass of 
80 GeV (see fig. 3, upper- right corner). A similar comparison for other model-cases also 
indicates that the constraints we obtain here are comparable to those one would derive from 
Fermi-LAT observations of local dwarf spheroidal galaxies {3^] • 
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FIG. 2: Left: The pair-annihilation cross section of singlet scalar dark matter into two photons. 
Solid lines correspond to a Lagrangian with only S 2 \H\ 2 interactions (62 = 0), and no relic abun- 
dance constraints, and a variety of values for the SU(2) Higgs mass nih = 120, 150, 180 GeV. 
Dashed lines correspond to a Lagrangian with S 2 \H\ 2 interactions plus a mass term for S (62 7^ 0), 
with specified fixed values of and 02, and again no relic abundance constraints. The dashed- 
dotted line features a quartic coupling, 02, fixed to satisfy relic abundance constraint. The dots 



correspond to the limits from the Fermi-LAT collaboration [17[ for different dark matter profiles, 
whereas the crosses are an extrapolation of such limits to lower energies. Right: constraints from 
Fermi data on the plane defined by rrih and A = (ms — m/ l /2)/m^. Here we set 62 — 0. 



200 



9 




100 120 140 160 180 20 



(a) 




20 40 60 80 100 120 140 160 180 200 



m s [GeV] 




1 1 1 1 1 1 1 1 1 1 1 1 mill 

20 40 60 80 100 120 140 160 180 200 
m [GeV] 



(b) 



o 




80 100 120 140 160 180 200 



(c) (d) 

FIG. 3: Top: Regions, on the (ms,mi,) parameter space, favored by the two events above back- 
ground observed by CDMS at 78% C.L. (green), and the regions excluded, at 90% C.L., also by 
CDMS (grey); we indicate the region ruled out by LEP (orange), favored by the S relic density 
(cyan), and excluded by Fermi searches for the monochromatic annihilation line, here coming from 
the process SS — > 77 (b 



limits from XENON 



38 



ack). In the left panel we set 62 = and we also show the controversial 
39I ]. while in the right panel we have 62 7^ and fix 02 = 0.1. Bottom: As 
in Fig. |3(b)[ but for a 2 = 0.05 (left) and 0.02 (right). 
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IV. IMPACT ON DIRECT DETECTION RESULTS 



In this section we explore how constraints on the singlet model from the annihilation 
SS 77 compare with the direct detection constraints. Only three parameters are relevant 
to our analysis: the Higgs mass m^, b 2 (or, alternatively, mg) and (Z2 1 . The direct detection 
constraints come from considering the spin-independent S'-nucleon cross section 

a\m%f 2 



^si 



2 4 ' 



(11) 



where mjy is the nucleon mass, / is the form factor. For numerical purposes, we set / = 1/3, 



following here Ref. 



21 



40). 



We will carry out the analysis for two different cases: (i) b 2 = and (ii) b 2 ^ 0. 



In this case, from Eq. ([3]) we have a 2 = m 2 s /v 2 , and we are left with only two free 
parameters, that can be traded off for the two particle masses and 171$. The regions 



excluded by LEP [411 ]. by CDMS [42j and by the Fermi results are shown in Fig. 3(a) . In the 
same plot, we also show the region compatible with the relic abundance and the one favored 



by CDMS at 78% confidence level 2l|, |42j . Notice that while the excluded region is obtained 



from the 90% C.L. upper limits on the WIMP-nucleon spin-independent cross section, the 
favored region comes from the two events observed by CDMS and is obtained at a different 
C.L. (namely, 78%). We see that there is some overlap between the two regions, which 
should not be regarded as an inconsistency in our approach. Such an overlap would in fact 
change if we used different choices for the confident levels. We also indicate the controversial 



limits from XENON 



3s 



391 ] with a dashed brown line. We notice that the monochromatic 



photon line limits are competitive with respect to the direct detection limits for large Higgs 
masses and for rris — rrih/2. We also find a small portion of the parameter space compatible 
with the tentative signal observed by CDMS that is ruled out by the monochromatic photon 
lines limit. The line constraints, however, never overlap with the region where S is thermally 
produced with the right relic abundance. 



1 Note that the singlet self quartic coupling 64 is completely irrelevant here. 
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FIG. 4: Regions excluded (at 90% C.L., brown) and regions marginally favored (at 78% C.L., 
green) by the CDMS results and regions excluded by Fermi constraints (black) on the SS — > 77 
process on the (7715,02) parameter space. Within the cyan region the S thermal relic abundance 
is compatible with the observed cosmological dark matter density; we also indicate the curves 
corresponding to a relic abundance of 0.01 (orange dashed lines) and of 0.001 (magenta dot-dashed 
lines). In the left panel we set rrih = 120 GeV, while in the right panel rrih = 200 GeV. 



B. b 2 + 0. 



Without restrictions on b 2 we have to deal with a three-parameter space: rrih, rns, 0,2- We 
first consider the plane (ms, rrih) and we show the excluded regions, as well as the CDMS 



favored region, for three different values of a 2 in Fig. 3(b) (where we set a 2 = 0.1), 3(c 



(a 2 = 0.05) and [3(d)] (a 2 = 0.02). 

Decreasing the value of a 2 squeezes the width of the resonance, as a smaller value for a 2 
reduces the cross section into two photons. At the same time, though, a reduced singlet-Higgs 
coupling also suppresses significantly the constraints from direct detection experiments. For 
a 2 < 0.05 and for ms > 60 GeV, the only constraints on the theory for viable values of the 
Higgs mass originate in fact from the two-photon annihilation mode. 

Lastly, we study the constraints on the plane (mj, a 2 ) in Fig. HJ where we fix rrih = 120 
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GeV (a) and rrih = 200 GeV (b). This cross-section of the theory parameter space illustrates 
that the SS —> 77 process excludes portions of the parameter space compatible with the 
CDMS putative signal (panel a), and it extends to values of smaller than those constrained 
by direct detection experiments, especially at growing values of the SM Higgs (panel b). 



We note that for > raw three body final states produced by annihilation into WW* 43 1 



as well as radiative electro-weak corrections 



44j can also play a significant role. 



V. THE ROLE OF THE VACUUM STABILITY CONSTRAINT 

The parameter space regions where the SS — > 77 annihilation mode puts further con- 
straints upon the scalar singlet dark matter model is broadened by examining the stability 
of the Higgs vacuum. The tree-level potential given by Eq. (T5]) can in fact easily develop a 
second minimum in the singlet direction in addition to the expected minimum in the Higgs 
direction. This happens when the mass-squared term 62 is negative and the coupling ci2 
is large enough to overcome the negative mass-squared of the Higgs field. Specifically, the 
singlet will have a minimum at S 2 = —62/^4 provided that a<i ■ — 2/i 2 > 0, which prevents 
the minimum from instead being a saddle point. If a second minimum exists and has a 
lower potential value than the minimum in the Higgs direction, then the physical vacuum 
state is unstable. However, there is a finite probability to tunnel to the new vacuum state. 
When the lifetime of the Higgs vacuum is longer than that of the universe, the vacuum is 



metastable and the theory is saved. We perform here a stability analysis similar to Ref. [301 ] . 
but extend their study to allow for metastability and the possibility for one-loop corrections 
to lead to instabilities in the effective potential. 

An important preliminary check is that the universe ends up in what will eventually 
evolve into the correct electro-weak minimum, at the electro-weak phase transition at high 
temperatures, instead of breaking the electro-weak symmetry in the singlet direction, and 
hence in what would then evolve to be the wrong zero-temperature vacuum. Finite temper- 
ature corrections to the effective potential tend to lift the potential away from the origin, 
and the lifting is strongest in regions with high mass particles. Even though we assume that 
the singlet scalar dark matter particle is relatively light in the electroweak phase, it could 
be quite heavy at the singlet minimum. In order to get a small dark matter mass either 
02 needs to be very small, in which case the electroweak minimum tends to be the true 



13 



vacuum (see Fig. [6] below), or -b 2 needs to be very large, which gives rise to large masses 
in the singlet direction. Therefore, it is reasonable to expect that at temperatures near 
the electroweak phase transition the finite temperature potential lifts the singlet minimum 
above the electroweak minimum, and the universe gets stuck in the soon to be metastable 
electroweak phase. A detailed analysis of the finite-temperature vacuum structure of the 
theory lies however beyond the scope of this analysis. 

The one loop corrections to the tree-level potential at zero temperature are 



=^±^™t(H°,S) [log (^) 



(12) 



where the sum is over all particle species, is the number of degrees of freedom per particle, 
rrii is the field-dependent particle mass, q = 3/2 for fermions and scalars and 5/6 for gauge 
bosons, A is the renormalization scale which we choose to be 1 TeV and H° = h + v indicates 
the neutral real component of the SU(2) complex doublet H. The (tree-level) Higgs mass is 
m 2 = /j> 2 + 3\(H°) 2 + a 2 S 2 , and the singlet mass is m| = b 2 + 3b 4 S 2 + a 2 (H ) 2 . Reference Q 
lists all other relevant field-dependent masses in its appendix. Note that they use slightly 
different notation than we do. They use A/6, m 2 and h where we use A, fi 2 , and H°. We 
follow the same procedure that they use for calculating the physical masses Mh and Ms- 

In order to ensure that the tree level minimum remains a minimum, we need to find the 
potential's second derivative: 

d 2 V x m 



OS 2 



2.32^^) [log ^ A2 ) 



1 



(13) 



Here, we took q = 3/2 and dropped terms containing d(m 2 )/dS which are zero at S = 
0. Unless the cutoff scale is taken to be smaller than the particle masses, the one-loop 
contribution tends to move the electroweak minimum towards instability. Usually this effect 
is not large enough to cancel the positive second derivative in the tree-level potential, but 
it can lead to instabilities in large sections of parameter space, as we will see below. 

To examine the problem of met ast ability, we must calculate the tunneling rate per unit 
volume from the metastable to stable vacua. This rate has the form F/V = Ae~ SE , where 



Se is the four- dimensional Euclidean action (see Ref. [45j for original work on the calculation 
of tunneling rates in field theory). The pref actor A is generally difficult to calculate, but its 
exact value matters little in comparison to the rate's reliance upon Se, so we can obtain an 
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approximate solution on dimensional grounds. Assuming an 0(4) symmetry in the equations 
of motion, the action is 



S r = 2vr 2 / r 3 dr 



(14) 



where r is the Euclidean coordinate r = (x,it). Minimizing Se produces the equations of 
motion 

and similarly for H° — > S, with the boundary conditions 

dH° dS . , 

— — = — = at r = 0; 16 
dr dr 

(H°, S) = (v, 0) at r = oo. (17) 

These equations describe a bubble of stable vacuum at r = embedded in a sea of the 
metastable electroweak vacuum. 

If there were only one field, then Eq. (Tl5l) could easily be solved by the under- 
shoot/overshoot method. One can exchange radius for time and then imagine the equation 
as describing a particle moving in the inverted potential — V(H°) with a peculiar 'time'- 
dependent friction term -(^j- + *fn). The particle starts near the absolute maximum of —V 
corresponding to the true vacuum, rolls down the potential, and then goes up again towards 
the maximum corresponding to the false vacuum. If the particle goes past the false vacuum, 
the initial conditions overshot the final conditions and they must be adjusted downward on 
the inverted potential. Conversely, if it does not make it to the false vacuum, the initial 
conditions undershot the final conditions and must be adjusted upwards. 

The two dimensional case, however, is much more complicated. We can simplify it by 
assuming that tunneling occurs along a fixed path parametrized by its path length: H° = 
H°(x), S = S(x), and (dH°/dx) 2 + (dS/dx) 2 = 1. Equation [151 then simplifies to the 
one-dimensional case 

^^ = >W,* (18) 
which we can solve by the undershoot /overshoot method. The trick then, is to choose 
the correct path. We do this by introducing a novel method of path deformation (see e.g. 



Rcf. 



46[ for another approach to finding the action. We will provide greater detail of our 



numerical algorithm in an upcoming paper.). 
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FIG. 5: Deformation of paths to solve the Euclidean equations of motion for Mh = 120 GeV, 
M s = 100 GeV, and o 2 = b A = 0.4. The electroweak minimum is at (H°,S) = (246,0) GeV, 
and the stable singlet minimum is (H°,S) = (0,197) GeV. Left: We choose both a straight line 
and an elliptical path as first guesses for the direction of tunneling. Arrows denote the 'normal' 
forces along each path. Right: The two guesses converge towards the correct solution. Each line 
represents 20 deformations with fixed step size. 

Let <p = (H°, S) describe the field coordinates. We can break the equations of motion 
into directions parallel and perpendicular to the direction of motion: 

f{^ + ~)=V*V\ (19) 
ax \dr z r dr J 11 

— * 

Then, imagining 0(x(r)) as describing a particle moving on a fixed track, the quantity 



(21) 



dx 2 \dr J 

corresponds to the normal force exerted by the track upon the particle. If the track coincides 
with the natural direction of motion, the normal force will be zero. Otherwise the normal 
force will point in the direction of necessary path deformation (see Fig. [5]). 

To execute the deformation, we first solve the one- dimensional equation of motion along 
a straight line between the two minima. We find dx/dr at n po i nt s = 100 evenly spaced points 
along the path, and then use this to find the normal force at those points. Each point 
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deforms an amount A<j) = aLN/\ W\ max , where a = 0.002 is our effective step size, L is the 
length of the path, and | Vy| ma x is the maximum absolute gradient of the potential along the 
path. A more rigorous approach would be to use an adaptive step size, but a small constant 
step size is sufficient for our purposes. Typically, the deformation converges onto a new path 
in roughly 100 steps, at which point we re-solve the one-dimensional equation of motion. 
In our cases, we only need to repeat this process two or three times before we achieve an 
accuracy of about 1% in the value of the Euclidean action. We check for convergence by 
picking a second starting path that lies on the other side of the final path (for example, an 
elliptical arc that connects the minima in the singlet and Higgs directions) and deforming 
from that direction. 

All that is left is to approximate the pre-factor A and find the critical value of Se for 
which we would have expected to see a phase transition. Here, we follow the argument in 

n 

Ref. [47J. Working in units of the electroweak scale, we set A = 1. The lifetime of the 
universe in electroweak units is e 101 , and the fraction of the universe filled with stable phase 
as a function of time is 1 — exp(— -j^t 4 ) (see Ref. 48[). Therefore, in order for the Higgs 
vacuum to be metastable we require that the action Se be greater than 404. 



We present our results in Fig. |6j Ignoring the one 



oop unstable region (red), we find 
identical areas of absolute stability as those in Ref. 30]. However, the total viable region 
of parameter space expands substantially by considering long-lived metastable vacua (blue 
regions). The metastable states add roughly 0.1 and 0.2 to the maximum allowed 02 value 
for low and high mass singlets, nearly doubling the allowed parameter space for theories 
with small singlet self-couplings (64). In the Figure, we show the regions compatible with 
a thermal relic abundance with a grey band, and points where the dark matter abundance 
equals 0.01 and 0.001 with dashed orange and dot-dashed magenta lines, respectively. We 
notice, by inspecting Fig. |4] and comparing with Fig. |6] that specifically for small singlet 
self-coupling 64, the region compatible with the putative CDMS signal interestingly overlaps 
largely with what we find are meta-stable but long-lived electro-weak vacuum configurations, 
for both heavy and light electro-weak masses. 

We super-impose the constraints we obtain from Fermi observations and the calculation 
of the cross section for the annihilation process SS — > 77 we presented here (black hatched 
areas). We notice that, especially for low values of 64, the two monochromatic photons anni- 
hilation mode sets strong constraints on the portion of the theory parameter space where the 
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-b, = 0.050 / Kg 
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FIG. 6: Regions of stability for different parameters in the scalar singlet dark matter model. White 
regions are absolutely stable; that is, the minimum at (H°,S) = (v,0) is the true vacuum. Blue 
regions are metastable with lifetimes longer than that of the universe, while green regions are 
metastable with lifetimes shorter than that of the universe. In red regions, the field configuration 
corresponding to (v, 0) is not a minimum. We show the regions compatible with a thermal relic 
abundance with a grey band, and points where the dark matter abundance equals 0.01 and 0.001 
with dashed orange and dot-dashed magenta lines, respectively. Finally, the black hatched areas 
indicate regions excluded by the SS — > 77 constraint. 



electro-weak vacuum is meta-stable with very long lifetimes. These regions also overlap with 
those compatible with the CDMS signal. Part of the meta-stable vacuum parameter space 
is nonetheless still open and not constrained by gamma-ray or direct detection experiments. 



VI. CONCLUSIONS 



We calculated the pair-annihilation cross section of a real, ^-symmetric scalar singlet ex- 
tension of the SM into two photons, and we derived the constraints on the theory parameter 
space from the Fermi limits on monochromatic gamma-ray lines. We showed that this new 
class of constraints improve on limits from direct dark matter searches in certain regions of 
parameter space, especially where the singlet mass is close to half the SM Higgs mass. The 
limits we find rule out portions of the theory parameter space compatible with the tentative 
signal events found by the CDMS Collaboration. We also calculated the lifetime of unstable 
vacuum configurations in the scalar potential, and showed that the gamma-ray limits are 
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quite relevant in regions where the electro- weak vacuum is meta-stable with a lifetime longer 
than the age of the universe. Those same regions also overlap with the portion of the theory 
parameter space compatible with the putative CDMS direct dark matter detection signal. 
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